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in this survey is a discretization of curvature line parametrized surfaces in Lie 
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cretization of curvature line parametrization which unifies the circular and conical 
nets. 
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1 Introduction 



The new field of discrete difi'erential geometry is presently emerging on the 
border between differential and discrete geometry. Whereas classical differ- 
ential geometry investigates smooth geometric shapes, discrete differential 
geometry studies geometric shapes with finite numbers of elements and aims 
to develop discrete equivalents of the geometric notions and methods of clas- 
sical differential geometry. The latter appears then as a limit of refinements 
of the discretization. Current interest in this field derives not only from 
its importance in pure mathematics but also from its relevance for com- 
puter graphics. An important example one should keep in mind here are 
polyhedral surfaces approximating smooth surfaces. 

One may suggest many different reasonable discretizations with the same 
smooth limit. Which one is the best? From the theoretical point of view 
the best discretization is the one which preserves all fundamental proper- 
ties of the smooth theory. Often such a discretization clarifies the struc- 
tures of the smooth theory and possesses important connections to other 
fields of mathematics (projective geometry, integrable systems, algebraic 
geometry, complex analysis etc.). On the other hand, for applications the 
crucial point is the approximation: the best discretization is supposed to 
possess distinguished convergence properties and should represent a smooth 
shape by a discrete shape with just few elements. Although these theo- 
retical and applied criteria for the best discretization are completely differ- 
ent, in many cases natural "theoretical" discretizations turn out to possess 
remarkable approximation properties and are very useful for applications 
[B^ILFWYW| . 

This interaction of the discrete and smooth versions of the theory led 
to important results in the surface theory as well as in the geometry of 
polyhedra. Classical achievements of discrete differential geometry are the 
fundamental results of Alexandrov and Pogorelov on metric geometry of 
polyhedra and convex surfaces: Alexandrov's theorem |Alj states that any 
abstract convex polyhedral metric is uniquely realized by a convex polyhe- 
dron in Euclidean 3-space. Pogorelov proved Pog the corresponding exis- 
tence and uniqueness result for abstract convex metrics by approximating 
smooth surfaces by polyhedra. 

Simplicial surfaces, i.e., discrete surfaces made from triangles, are ba- 
sic in computer graphics. This class of discrete surfaces, however, is too 
unstructured for analytical investigation. An important tool in the theory 
of smooth surfaces is the introduction of (special) parametrizations of a 
surface. Natural analogues of parametrized surfaces are quadrilateral sur- 
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faces, i.e. discrete surfaces made from (not necessarily planar) quadrilat- 
erals. The strips of quadrilaterals obtained by gluing quadrilaterals along 
opposite edges are analogs of coordinate lines. Probably the first nontrivial 
example of quadrilateral surfaces studied this way are discrete surfaces with 
constant negative Gaussian curvature introduced by Sauer and Wunderlich 
|Sll IWuj . Currently discrete parametrized surfaces are becoming more im- 
portant in computer graphics. They lead to meshes that better represent 
the shape of the surface and look regular jACSPLDl IDKHI IMKI iLPWYWj . 

It is well known that differential equations describing interesting spe- 
cial classes of surfaces and parametrizations are integrable (in the sense of 
the theory of integrable systems), and, conversely, many of interesting inte- 
grable systems admit a differential-geometric interpretation. A progress in 
understanding of the unifying fundamental structure the classical differential 
geometers were looking for, and simultaneously in understanding of the very 
nature of integrability, came from the efforts to discretize these theories. It 
turns out that many sophisticated properties of differential-geometric ob- 
jects find their simple explanation within the discrete differential geometry. 
The early period of this development is documented in the work of Sauer |S2j . 
The modern period began with the work by Bobenko and Pinkall |BPlllBP2] 
and by Doliwa and Santini |DS|I(]D^ . A closely related development of the 
spectral theory of difference operators on graphs was initiated by Novikov 
with collaborators |ND| INH IN2j , see also |DN : for a further development of 
a discrete complex analysis on simplicial manifolds. 

Discrete surfaces in Euclidean 3-space is the basic example considered in 
this survey. This case has all essential features of the theory in all generality, 
generalizations for higher dimensions are straightforward. On the other 
hand, our geometric three-dimensional intuition helps to understand their 
properties. 

Discrete differential geometry related to integrable systems deals with 
multidimensional discrete nets, i.e., maps from the regular cubic lattice Z*" 
into M'^ specified by certain geometric properties (as mentioned above, we 
will be most interested in the case = 3 in this survey) . In this setting dis- 
crete surfaces appear as two dimensional layers of multidimensional discrete 
nets, and their transformations correspond to shifts in the transversal lattice 
directions. A characteristic feature of the theory is that all lattice direc- 
tions are on equal footing with respect to the defining geometric properties. 
Discrete surfaces and their transformations become indistinguishable. We 
associate such a situation with the multidimensional consistency, and this is 
one of our fundamental discretization principles. The multidimensional con- 
sistency, and therefore the existence and construction of multidimensional 
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nets, relies just on certain incidence theorems of elementary geometry. 

Conceptually one can think of passing to a continuum limit by refin- 
ing mesh size in some of the lattice directions. In these directions the net 
converges to smooth surfaces whereas those directions that remain discrete 
correspond to transformations of the surfaces (see Fig. P). The smooth 



Figure 1: From the discrete master theory to the classical theory: surfaces 
and their transformations appear by refining two of three net directions. 

theory comes as a corollary of a more fundamental discrete master theory. 
The true roots of the classical surface theory are found, quite unexpectedly, 
in various incidence theorems of elementary geometry. This phenomenon, 
which has been shown for many classes of surfaces and coordinate systems 
|BP2| IBS2j , currently is getting accepted as one of the fundamental features 
of classical integrable differential geometry. 

Note that finding simple discrete explanations for complicated differen- 
tial geometric theories is not the only outcome of this development. Having 
identified the roots of the integrable differential geometry in the multidi- 
mensional consistency of discrete nets, we are led to a new (geometric) 
understanding of the integrability itself |BS1L lABSL IBS2| . 

The simplest and at the same time the basic example of consistent mul- 
tidimensional nets are multidimensional Q-nets |DSj . or discrete conjugate 
nets |S2j . which are characterized by planarity of all quadrilaterals. The pla- 
narity property is preserved by projective transformations and thus Q-nets 
are subject of projective geometry (like conjugate nets, which are smooth 
counterparts of Q-nets). 

Here we come to the next basic discretization principle. According to 
F. Klein's Erlangen program, geometries are classified by their transforma- 
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tion groups. Classical examples are projective, affine, Euclidean, spherical, 
hyperbolic geometry, and sphere geometries of Lie, Moebius, and Laguerre. 
We postulate that the transformation group as the most fundamental fea- 
ture should be preserved by a discretization. This can be seen as a sort of 
discrete Erlangen program. 

Thus we come to the following fundamental 

Discretization Principles: 

• Transformation group principle: smooth geometric objects and their 
discretizations belong to the same geometry, i.e. are invariant with 
respect to the same transformation group. 

• Consistency principle: discretizations of smooth parametrized geome- 
tries can be extended to multidimensional consistent nets. 

Let us explain why such different imperatives as the transformation 
group principle and the consistency principle can be simultaneously im- 
posed for discretization of classical geometries. The transformation groups 
of various geometries, including those of Lie, Mobius and Laguerre, are sub- 
groups of the projective transformation group. Classically, such a subgroup 
is described as consisting of projective transformations which preserve some 
distinguished quadric called absolute. A remarkable result by Doliwa |Dlj 
is that multidimensional Q-nets can be restricted to an arbitrary quadric. 
This is the reason why the Discretization Principles work for the classical 
geometries. 

In this survey we deal with three classical geometries described in terms 
of spheres: Mobius, Laguerre and Lie geometries. They have been developed 
by classics, the most elaborated presentation of these geometries can be 
found in Blaschke's book IBl_. 

Mobius geometry is the most popular one of these three geometries. It 
describes properties invariant with respect to Mobius transformations which 
are compositions of reflections in spheres. For > 3, Mobius transforma- 
tions of coincide with conformal transformations. Mobius geometry does 
not distinguish between spheres and planes (planes are regarded as spheres 
through the infinitely remote point oo, which compactifies to the N- 
sphere S^). On the other hand, points are considered as objects different 
from spheres. Surfaces are described through their points. Classical exam- 
ples of Mobius-geometric properties of surfaces are conformal parametriza- 
tion and the Willmore functional |Wij . Recent progress in this field is to 
a large extent due to interrelations with the theory of integrable systems 
[FLPHlTj. 
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Laguerre geometry does not distinguish points and spheres (points are 
treated as spheres of zero radius). On the other hand, planes are distin- 
guished. Surfaces are described through their tangent planes. A particular 
Laguerre transformation of a surface is a shift of all tangent planes in the 
normal direction at a constant distance. This transformation is called the 
normal shift. 

Lie geometry is a natural unification of Mobius and Laguerre geometries: 
points, planes and spheres are treated on an equal footing. The transfor- 
mation group is generated by Mobius transformations and the normal shift 
transformations. Surfaces are described through their contact elements. A 
contact element can be understood as a surface point together with the cor- 
responding tangent plane. The one-parameter family of spheres through a 
point with a common tangent plane gives a Lie-geometric description of a 
contact element. The point of the surface and the tangent plane at this 
point are just two elements of this family. 

Integrability aspects of the surface theory in Lie geometry have been 
studied by Ferapontov |Fll IF2j , Musso and Nicolodi |MN| , and Burstall and 
Hertrich-Jeromin |BH.TH lBH,T2j . 

The main concrete geometric problem discussed in this survey is a dis- 
cretization of curvature line parametrized surfaces. Curvature lines are in- 
tegral curves of the principal directions. Any surface away from its umbilic 
points can be parametrized by curvature lines. Curvature line parametriza- 
tion is attracting the attention of mathematicians and physicists for two 
centuries. The classical results in this field can be found in the books by 
Darboux |Dall IDa2j and Bianchi [Bij. In particular a classical result of 
Dupin [Dj claims that the coordinate surfaces of triply orthogonal coor- 
dinate systems intersect along their common curvature lines. Ribaucour 
has discovered a transformation of surfaces preserving the curvature line 
parametrization (see |E2,). A surface and its Ribaucour transform envelope 
a special sphere congruence. Bianchi has shown |Bi2j that Ribaucour trans- 
formations are permutable: given two Ribaucour transforms of a surface 
there exists a one-parameter family of their common Ribaucour transforms. 

Recently curvature line parametrizations and orthogonal systems came 
back into the focus of interest in mathematical physics as an example of an 
integrable system. Zakharov [Z] has constructed a variety of explicit solu- 
tions with the help of the dressing method. Algebro-geometric orthogonal 
coordinate systems were constructed by Krichever |Kj. The recent interest 
to this problem is in particular motivated by applications to the theory of 
the associativity equations developed by Dubrovin |Duj . Remarkable geo- 
metric properties make curvature line parametrizations especially useful for 
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visualization of surfaces in computer graphics |ACSDLDl ILPWYWj . 

The question of proper discretization of the curvature hne parametrized 
surfaces and orthogonal systems became recently a subject of intensive 
study. Circular nets, which are Q-nets with circular quadrilaterals, as dis- 
crete analogs of curvature line parametrized surfaces were mentioned by 
Nutbourne and Martin |NMj . Special circular nets as discrete isothermic 
surfaces were investigated in |BP2j . The circular discretization of triply- 
orthogonal coordinate systems was first suggested by one of the authors in 
[B] , Doliwa and Santini |DSj made the next crucial step in the development 
of the theory. They considered discrete orthogonal systems as a reduction of 
discrete conjugated systems |(yDSj . generalized them to arbitrary dimension 
and proved their multidimensional consistency based on the classical Miquel 
theorem |Bej . 

Matthes and the authors of this survey have proven |BMSj that circular 
nets approximate smooth curvature line parametrized surfaces and orthog- 
onal systems with all derivatives. Numerical experiments show that circular 
nets have the desired geometrical properties already at the coarse level and 
not only in the refinement limit as it approaches a smooth curvature line pa- 
rameterized surface. This is important for applications in computer graphics 
II^PWYWj. 

A convenient analytic description of circular nets has been given by 
Konopelchenko and Schief |KSj . Analytic methods of the soliton theory 
have been applied to circular nets by Doliwa, Manakov and Santini |DMSj 
(9-method) and by Akhmetshin, Volvovskii and Krichever |AKVj (algebro- 
geometric solutions). Bobenko and Hertrich-Jeromin IBHJ' have given a 
Clifford algebra description of circular nets. 

Circular nets are preserved by Mobius transformations, and thus should 
be treated as the discretization of curvature line parametrizations in Mobius 
geometry. A remarkable recent development by Liu, Pottmann, Wallner, 
Yang, and Wang |LPWYW| is the introduction of conical nets, which should 
be treated as the discretization of curvature line parametrizations in La- 
guerre geometry. These are special Q-nets characterized by the property 
that four quadrilaterals meeting at a vertex are tangent to a common cone 
of revolution. Equivalently, conical nets can be characterized as Q-nets with 
circular Gauss maps, i.e., the unit normals to the quadrilaterals comprise a 
circular net in the unit sphere S'^. Circular Gauss maps defined at vertices 
of a given circular net were previously introduced by Schief |Sch2j . however 
without relation to conical nets. Conical nets, like circular ones, satisfy the 
second discretization principle (consistency). 

In the present survey, we find a discretization of curvature line parametriza- 
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tion which unifies the circular and the conical nets by systematically apply- 
ing the Discretization Principles. 




Figure 2: Principal directions through touching spheres. 

It is well known that curvature lines are subject of Lie geometry, i.e., are 
invariant with respect to Mobius transformations and normal shifts. To see 
this, consider an infinitesimal neighborhood f/ of a point x of an oriented 
smooth surface in M^, and the pencil of spheres S{r) of the signed radii r, 
touching the surface at x, see Fig.[2j The signed radius r is assumed positive 
if S{r) lies on the same side of the surface as the normal n, and negative 
otherwise; 5(oo) is the tangent plane. For small ro > the spheres S'(ro) 
and S'(— ro) intersect U in x only. The set of the touching spheres with 
this property (intersecting U in x only) has two connected components: 
containing S{ro) and M__ containing rg) for small rg > 0. The boundary 
values 

ri = sup{r : S{r) € M+j, r2 = inf{r : S{r) € M_} 

are the principal curvatures of the surface in x. The directions in which 
S{ri) and S{r2) touch U are the principal directions. 

Clearly, all ingredients of this description are Mobius-invariant. Under a 
normal shift by the distance d the centers of the principal curvature spheres 
are preserved and their radii are shifted by d. This implies that the principal 
directions and thus the curvature lines are preserved under normal shifts, as 
well. 

A Lie-geometric nature of the curvature line parametrization yields that 
it has a Lie-invariant description. Such a description can be found in 
Blaschke's book |Blj . A surface in Lie geometry, as already said, is con- 
sidered as consisting of contact elements. Two infinitesimally close contact 
elements (sphere pencils) belong to the same curvature line, if and only if 
they have a sphere in common, which is the principal curvature sphere. 
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By a literal discretization of this Blaschke's Lie-geometric description of 
smooth curvature line parametrized surfaces, we define a discrete principal 
contact element net as a map I? {contact elements of surfaces in M^} 
such that any two neighboring contact elements have a sphere in common. 

In the projective model of Lie geometry spheres in (including points 
and planes) are represented by elements of the so called Lie quadric L C MP^, 
contact elements are represented by isotropic lines, i.e., lines in L, sur- 
faces are represented by congruences of isotropic lines. In the curvature 
line parametrization, the parametric families of isotropic lines comprise de- 
velopable surfaces in L. 

Accordingly, a discrete principal contact element net in the projective 
model of Lie geometry is a discrete congruence of isotropic lines 

(. -.1?' ^ {isotropic lines in L} 

such that any two neighboring lines intersect. Intersection points of neigh- 
boring lines are, as in the smooth case, the principal curvature spheres. 
They are associated with the edges of Z^. Four principal curvature spheres 
associated to the edges with a common vertex belong to the same contact 
element, i.e., have a common touching point. 




Figure 3: Geometry of principal contact element nets. Four neighboring 
contact elements are represented by points and (tangent) planes. The points 
are concircular, the planes are tangent to a cone of revolution. Neighboring 
normal lines intersect at the centers of principal curvature spheres. 
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In projective geometry, discrete line congruences have been introduced 
by Doliwa, Mafias and Santini |DMSj . Discrete line congruences are closely 
related to Q-nets, and, like the latter, are multidimensionally consistent. It 
follows from our results that they can be restricted to the Lie quadric (actu- 
ally, to any ruled quadric). Thus, principal contact element nets satisfy the 
second discretization principle. In particular, this yields discrete Ribaucour 
transformations between principal contact element nets. 

The Lie-geometric notion of discrete principal contact element nets uni- 
fies the Mobius-geometric one (circular nets) and the Laguerre-geometric 
one (conical nets). Indeed, any contact element I contains a point x and a 
plane P . It turns out that for a surface 

i ■.Tl? ^ {isotropic lines in L} = {contact elements in M^}, 

the points comprise a circular net 

whereas the planes comprise a conical net 

P -.1? ^ {planes in M^}. 

The corresponding geometry is depicted on Fig. 13 Schematically, this Lie- 
geometric merging of the Mobius- and Laguerre-geometric notions is pre- 
sented on Fig. ^ 

This survey is organized as follows. In Sect. [2 we start with a review of 
the basic multidimensionally consistent systems - the Q-nets and discrete 
line congruences. The basic notions of Lie, Mobius and Laguerre geometries 
are briefly presented in Sect. 01 Sect. |1] contains main new results on discrete 
curvature line parametrized surfaces: Lie-geometric Definitions ^1 El and 
Theorem which describes interrelations of discrete curvature line nets in 
Lie, Mobius and Laguerre geometries. Geometric characterization of Ribau- 
cour transformations and discrete R-congruences of spheres as quadrilateral 
nets in the Lie quadric is given in Sect. El 

Let us note that, due to the classical Lie's sphere-line correspondence, 
the Lie-geometric theory presented in this survey can be transferred to the 
context of projective line geometry in three-space: the Lie quadric is re- 
placed by the Pliicker quadric, the curvature lines and R-congruences of 
spheres correspond to the asymptotic lines and the W-congruences of lines, 
respectively. The projective theory of discrete asymptotic nets has been 
developed by Doliwa JD2j. 
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Figure 4: Geometry of principal contact element nets. Four neighboring 
contact elements produce a hexahedron with vertices in the Lie quadric 
L and with planar faces. The bottom quadrilateral is the intersection of 
the three-dimensional space V = span(^, ^i, ^12) with the 4-space in MP^ 
representing points in M^. The top quadrilateral is the intersection of V 
with the 4-space in MP^ representing planes in M^. Each side quadrilateral 
lies in the plane of two intersecting lines ^ C L. 



Our research in discrete differential Lie geometry has been stimulated 
by the recent introduction of conical nets by Liu, Pottmann, Wallner, Yang, 
and Wang |LPWYWj . The advent of the second (after circular nets) dis- 
cretization of curvature line parametrizations posed a question about the 
relation between the different discretizations. Independently, a relation be- 
tween circular and conical nets has been found by Pottmann P . We are 
grateful to H. Pottmann and J. Wallner for numerous communications on 
conical nets and for providing us with their unpublished results. We thank 
also U. Pinkall for useful discussions. 

2 Consistency as a discretization principle 
2.1 Q-nets 

We use the following standard notation: for a function / on Z™" we write 

Tif{u) = f{u + ei), 

where Cj is the unit vector of the i-th coordinate direction, 1 < i < m. We 
use also the shortcut notations fi for Tif, fij for TiTjf, etc. 
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The most general of the known discrete 3D systems possessing the prop- 
erty of 4D consistency are nets consisting of planar quadrilaterals, or Q- 
nets. Two-dimensional Q-nets were introduced by Sauer |S2j . the multi- 
dimensional generalization has been given by Doliwa and Santini |DSj . Our 
presentation in this section follows the latter paper. The fundamental im- 
portance of multi-dimensional consistency of discrete systems as their inte- 
grability has been put forward by the authors |BSH IABS| IBS2| . 

Definition 1 (Q-net) A map f : — > RP^ is called an m- dimensional 
Q-net (quadrilateral net, or discrete conjugate net) in MP'^ [N > 3), if all 
its elementary quadrilaterals if,fi,fij,fj) (at any u € and for all pairs 
l^^T^i^^n-j o,re planar. 

Thus, for any elementary quadrilateral, any representatives /, /j, fij 
of its vertices in the space M^^^ of homogeneous coordinates satisfy an 
equation of the type 



Representatives in any hyperplane of M^+^, for instance, in the affine part 
of MP^ = P(]R^+i), satisfy such an equation with 1 = dj + Cji + pij, 



Given three points /, /i, /2 in MP , one can take any point of the plane 
through these three points as the fourth vertex /12 of an elementary quadri- 
lateral (/, /i, /i25 72) of a Q-net. Correspondingly, given any two discrete 
curves / : Z X {0} ^ MP^ and / : {0} x Z ^ MP^ with a common point 
/(0, 0), one can construct infinitely many Q-surfaces / : Z^ — > MP^ with 
these curves as coordinate ones: the construction goes inductively, on each 
step one has a freedom of choosing a point in a plane (two real parameters) . 

On the other hand, constructing elementary hexahedra of Q-nets cor- 
responding to elementary 3D cubes of the lattice Z™ admits a well-posed 
initial value problem with a unique solution, therefore one says that Q-nets 
are described by a discrete 3D system: 

Theorem 2 (Elementary hexahedron of a Q-net) Given seven points 
f, fi and fij (1 < « < J < 3) in MP^, such that each of the three quadrilat- 
erals (/, fi, fij, fj) is planar (i.e., fij lies in the plane Yiij through f , fi, fj), 
define three planes T^Jlij as those passing through the point triples fk, fik, 
fjk, respectively. Then these three planes intersect generically at one point: 




(1) 



that is 



fij - f = Cijifj - f) + Cji{fi - f). 



(2) 



/i23 = nn23 n T2ni3 n r3ni2 . 



12 



Proof. Planarity of the quadrilaterals {f,fi,fij,fj) assures that all seven 
initial points /, /, and fij belong to the three-dimensional space ni23 through 
the four points /, /i, /2, /s- Hence, the planes T^Hij lie in this three- 
dimensional space, and therefore generically they intersect at exactly one 
point. □ 



The elementary construction step from Theorem]^ is symbolically rep- 
resented on Fig. |SJ which is the picture we have in mind when thinking and 
speaking about discrete three-dimensional systems with dependent variables 
(fields) attached to the vertices of a regular cubic lattice. 

As follows from Theorem |21 a three-dimensional Q-net / : — > MP^ is 
completely determined by its three coordinate surfaces 

/ : X {0} ^ MP^, / : Z X {0} X Z ^ RF^ , / : {0} x ^ MP^. 

Turning to an elementary cube of the dimension m > 4, we see that one 
can prescribe all points /, fi and fij for all 1 < i < j < m. Indeed, these 
data are clearly independent, and one can construct all other vertices of an 
elementary cube starting from these data, provided one does not encounter 
contradictions. To see the possible source of contradictions, consider in 
detail first the case of m = 4. From /, fi and fij (1 < i < j < 4) one 
determines all fijk uniquely. After that, one has, in principle, four different 
ways to determine /1234, from four 3D cubic faces adjacent to this point; 
see Fig. |H1 Absence of contradictions means that these four values for /1234 
automatically coincide. We call this property the 4D consistency. 

Definition 3 (4D consistency) A 3D system is called 4D consistent, if it 
can he imposed on all three-dimensional faces of an elementary cube of Z^. 





Figure 5: 3D system on an elementary cube 
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Figure 6: 4D consistency of 3D systems 



Remarkably, construction of Q-ncts based on the planarity of all elemen- 
tary quadrilaterals enjoys this property. 

Theorem 4 (Q-nets are 4D consistent) The 3D system governing Q- 
nets is 4D- consistent. 

Proof. In the construction above, the four values in question are 

/l234 = TiT2n34 11X1X31124 1X1X41123 , 

and three other ones obtained by cyclic shifts of indices. Thus, we have to 
prove that the six planes XiXjll^.^ intersect in one point. 

First, assume that the ambient space MP^ has dimension > 4. Then, 
in general position, the space 111234 through the five points /, /i (1 < ^ < 4) 
is four-dimensional. It is easy to understand that the plane TiTjUke is the 
intersection of two three-dimensional subspaces TiUjki and TjUike- Indeed, 
the subspacc TiUjki through the four points /j, fij, fik, fu contains also fijk, 
fiji, and fjke- Therefore, both TiUjke and Tjllikt contain three points fij, 
fijk and fiji, which determine the plane TiTjIlk£. Now the intersection in 
question can be alternatively described as the intersection of the four three- 
dimensional subspaces xin234, X2ni34, X3ni24 and X4ni23 of one and the 
same four-dimensional space 111234. This intersection consists in the generic 
case of exactly one point. 

In the case of = 3, we embed the ambient space into MP^, then 
slightly perturb the point f^ by adding a small component in the fourth 
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coordinate direction, then apply the above argument, and after that send 
the perturbation to zero. This proof works since, as one can easily see, on 
each step of the construction the perturbation remains regular. □ 

The m-dimensional consistency of a 3D system for m > 4 is defined 
analogously to the m = 4 case. Remarkably and quite generally, the 4- 
dimensional consistency already implies m-dimensional consistency for all 
m > 4. 

Theorem 5 (4D consistency yields consistency in all higher dimen- 
sions) Any 4D consistent discrete 3D system is also m-dimensionally con- 
sistent for any m > 4. 

Proof goes by induction from the (m — l)-dimensional consistency to the m- 
dimensional consistency, but, for the sake of notational simplicity, we present 
the details for the case m = 5 only, the general case being absolutely similar. 

Initial data for a 3D system on the 5D cube C12345 with the fields on 
vertices consist of the fields /, fi and fij for all 1 < i < j < 5. From these 
data one first gets ten fields fijk foi 1 < i < j < k < 5, and then five fields 
fijke for l<i<j<k<i<5 (the fact that the latter are well defined 
is nothing but the assumed 4D consistency for the 4D cubes Qijki)- Now, 
one has ten possibly different values for /12345, coming from ten 3D cubes 
TiTj^kim- To prove that these ten values coincide, consider five 4D cubes 
TiQjkim- For instance, for the 4D cube T1C2345 the assumed consistency 
assures that the four values for /12345 coming from four 3D cubes 

''"IT2C345, nT3C245, nT"4C235, nT"5C234 

are all the same. Similarly, for the 4D cube r2Ci345 the 4D consistency leads 
to the conclusion that the four values for /12345 coming from 

n''"2C345, 7-2T3C145, 7"2T4Cl35, ''"2T5C134 

coincide. Note that the 3D cube tiT2Qm^, the intersection of riC2345 and 
^"261345) is present in both lists, so that we now have seven coinciding values 
for /12345. Adding similar conclusions for other 4D cubes TiQjkirm we arrive 
at the desired result. □ 

Theorems m 13 yield that Q-nets are m-dimensionally consistent for any 
m > 4. This fact, in turn, yields the existence of transformations of Q-nets 
with remarkable permutability properties. Referring for details to |DSMl 
IBS2) . we mention here only the definition. 
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Definition 6 (F-transformation of Q-nets) Two m- dimensional Q-nets 
f,f^ : Z*" MP^ are called F-transforms (fundamental transforms) of 
one another, if all quadrilaterals {f , fi, fi' , f^) (at any u G and for all 
1 < i < m) are planar, i.e., if the net F : x {0,1} MP^ defined 
by F(n, 0) = f{u) and F{u, 1) = f^{u) is a two-layer (m + 1)- dimensional 
Q-net. 

It follows from Theorem |2 that, given a Q-net /, its F-transform is 
uniquely defined as soon as its points along the coordinate axes are suitably 
prescribed. 

2.2 Discrete line congruences 

Another important geometrical objects described by a discrete 3D system 
which is 4D consistent, are discrete line congruences. Their theory has been 
developed by Doliwa, Santini and Mafias |DSM| . whose presentation we 
follow in this section. 

Let be the space of lines inRP^; it can be identified with the Grass- 
mannian Gr(A'^ + 1, 2) of two-dimensional vector subspaces of M^+^. 

Definition 7 (Discrete fine congruence) A map I : — > iL^ is called 
an m- dimensional discrete line congruence in MP^ [N > 3), if any two 
neighboring lines i, ii (at any u G and for any 1 < i < m) intersect (are 
co-planar). 

For instance, lines i = (ff^) connecting corresponding points of two 
Q-nets /, : Z™" —>■ WF^ in the relation of F-transformation clearly build 
a discrete line congruence. 

A discrete line congruence is called generic, if for any u £ Z™ and for 
any 1 < i ^ j ^ k ^ i < m, the four lines I, ii, ij and Ik span a four- 
dimensional space (i.e., a space of a maximal possible dimension). This 
yields, in particular, that for any u G Z™ and for any 1 < i ^ j < m, the 
three lines i, ii and ij span a three-dimensional space. 

Construction of line congruences is similar to that of Q-nets. Given 
three lines i, i\, i^ of a congruence, one has a two-parameter family of lines 
admissible as the fourth one i\2: connect by a line any point of i\ with any 
point of £2- Thus, given any two sequences of lines I : Ij x {0} — and 
i : {0} X Z — !■ such that any two neighboring lines are co-planar, one 
can extend them to a two-dimensional line congruence f : 1? ^ in an 
infinite number of ways: on each step of the inductive procedure one has a 
freedom of choosing a line from a two-parameter family. 
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The next theorem shows that non-degenerate hne congruences are de- 
scribed by a discrete 3D system: 

Theorem 8 (Elementary hexahedron of a discrete line congruence) 

Given seven lines £, £i and iij (1 < i < i < 3) in MP^, such that £ intersects 
each of ii, the space V123 spanned by i, ii, £2, is has dimension four, and 
each £i intersects both iij and in., there is a unique line ^123 that intersects 
all three iij. 

Proof. All seven lines, and therefore also the three-dimensional spaces 
TiVjk = span{ii,iij,iik) lie in V123. A line that intersects all three of iij 
should lie in the intersection of these three three-dimensional spaces. But a 
generic intersection of three three-dimensional spaces in V123 is a line: 

423 = riV23 n T2V13 n T3V12. 

It is now not difficult to realize that this line does, indeed, intersect all three 
of iij. For instance, T1V23 n r2Vi3 = span(^i2, £13) n span(^i2, ^23) is a plane 
containing ii2, therefore its intersection with T3V12 (the line £123) intersects 

il2. □ 

A similar argument shows: 

Theorem 9 (Discrete line congruences are 4D consistent) The 3D 

system governing discrete line congruences is 4^D- consistent. 

Like in the case of Q-nets, this theorem yields the existence of transfor- 
mations of discrete line congruences with remarkable permutability proper- 
ties. 

Definition 10 (F-transformation of line congruences) Two m-dimen- 
sional line congruences i, i'^ : Z™ are called F-transforms of one 

another, if the corresponding lines i and i^ intersect (at any u € 17^), i.e., if 
the map L : Z"" x {0, 1} defined by L{u, 0) = i{u) and L{u, 1) = i+{u) 

is a two-layer (m + 1)- dimensional line congruence. 

Again, it follows from Theorem [J] that, given a line congruence i, its F- 
transform i^ is uniquely defined as soon as its lines along the coordinate 
axes are suitably prescribed. 

According to Definitional any two neighboring lines i = t{u) and ii = 
i{u + ei) of a line congruence intersect at exactly one point f = inii G RP^ 
which is thus combinatorially associated with the edge (n, n + e^) of the 
lattice Z"*: / = /(n, n + Cj). It is, however, sometimes more convenient to 
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Figure 7: Four lines of a congruence 

use the notation f{u,u + Cj) = f^^\u) for this points, thus associating it to 
the vertex u of the lattice (and, of course, to the coordinate direction i). See 
Fig. El 

Definition 11 (Focal net) For a discrete line congruence £ : HT^ , 
the map : Z"^ ^ MP^ defined by f^\u) = £{u) D £{u + a) is called its 
i-th focal net. 

Theorem 12 For a non- degenerate discrete line congruence £ : Z*" — > , 
all its focal nets f^^^ : 11^ MP^, \<k<ra, are Q-nets. 

Proof consists of two steps. 

► First, one shows that for the k-th focal net f^''\ all elementary quadri- 
laterals {f^^\ ff'\ f-j^\ f'^'^) are planar. This is true for any line congruence. 
Indeed, both points /^'^^ and Z^^'* lie on the line £k, while both points f[^'^ and 

f^^ lie on the line £ik- Therefore, all four points lie in the plane spanned by 
these two lines £k and £ik which intersect by definition of a line congruence. 

► Second, one shows that for the fc-th focal net f^^\ all elementary 
quadrilaterals {f^^\ fl'^\ f[^\ ff^'^), with both i j different from A;, are 
planar. Here, one uses essentially the assumption that the line congruence 
£ is generic. All four points in question lie in each of the three-dimensional 
spaces 

Vij = span(^, £i, £j,£ij) and TkVij = span{£k, £ik,£jk,£ijk) 
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(see Fig.lHI). Both 3-spaces lie in the four-dimensional space Vijk = span(£, £i, £j 
so that generically their intersection is a plane. □ 




Figure 8: Elementary (ij) quadrilateral of the k-th focal net 

Corollary 13 (Focal net of F-transformation of a line congruence) 

Given two generic line congruences i, : in the relation of F- 

transformation, the intersection points f = i (1 form a Q-net f : HT^ — > 

2.3 Q-nets in quadrics 

We consider an important admissible reduction of Q-nets: they can be con- 
sistently restricted to an arbitrary quadric in MP^. In the smooth differen- 
tial geometry, i.e., for conjugate nets, this is due to Darboux |Dalj . In the 
discrete differential geometry this result has been found by Doliwa |Dlj . 

A deep reason for this result is the following fundamental fact well known 
in classical projective geometry (see, e.g., |B2l ): 

Theorem 14 (Associated point) For any seven points of CP'^ in gen- 
eral position, there exists the eighth point (called the associated one), which 
belongs to any quadric through the original seven points. 
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Proof is based on the following computations. The equation Q = of a 
quadric in CP^ has ten coefficients (homogeneous polynomial of 4 variables). 
Therefore, a unique quadric Q = can be drawn through nine points in gen- 
eral position. Similarly, a pencil (one-parameter linear family) of quadrics 
Q -|- AQ' = can be drawn through eight points in general position, and a 
two-parameter linear family of quadrics Q -|- AQ' -|- /iQ" = can be drawn 
through seven points in general position. Generically, solution of a system 
of three quadratic equations 

Q = 0, Q' = 0, Q" = 

for the intersection of three quadrics in CP^ consists of eight points. It 
can be shown that the three quadrics spanning the above-mentioned two- 
parameter family can be considered generic enough for such a conclusion. 
Clearly, the resulting eight points lie on every quadric of the two-parameter 
family. □ 

Theorem 15 (Elementary hexahedron of a Q-net in a quadric) If 

seven points f, fi, and fij (1 < i < j < 3) of an elementary hexahedron 
of a Q-net f -.11^ ^ MP^ belong to a quadric Q C MP^, then so does the 
eighth point /123 . 

Proof. The original seven points can be considered lying in a three-dimen- 
sional space, and they are known to belong to three (degenerate) quadrics - 
the pairs of planes Yljk^TiUjk for (jfc) = (12), (23), (31). Clearly, the eighth 
intersection point of these quadrics is /123 = Tin23 nr2n3i nr3ni2, and this 
has to be the associated point. According to Theorem [TH it belongs to any 
quadric through the original seven points, in particular, to Q. □ 

3 Geometries of spheres 
3.1 Lie geometry 

A classical source on Lie geometry is Blaschke's book |B1) . see also a modern 
account by Cecil [C]. 

Following geometric objects in the Euclidean space M.^ are elements of 
Lie geometry: 

• Oriented hyperspheres. A hypersphere in M.^ with center c G M.^ and 
radius r > is described by the equation S = {x £ : |x — cp = r^}. 
It divides in two parts, inner and outer. If one denotes one of 
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two parts of M as positive, one comes to the notion of an oriented 
hypersphere. Thus, there are two oriented hyperspheres for any 
S. One can take the orientation of a hypersphere into account by 
assigning a signed radius itr to it. For instance, one can assign positive 
radii r > to hyperspheres with the inward field of unit normals and 
negative radii r < to hyperspheres with the outward field of unit 
normals. 

• Oriented hyperplanes. A hyperplane in M.'^ is given by the equation 
P = {x £ : {v,x) = d}, with a unit normal v G S^~-^ and 
d € M. Clearly, the pairs {v,d) and {—v,—d) represent one and the 
same hyperplane. It divides in two halfspaces. Denoting one of 
two halfspaces as positive, one arrives at the notion of an oriented 
hyperplane. Thus, there are two oriented hyperplanes for any P. 
One can take the orientation of a hypersphere into account by assigning 
the pair (f , d) to the hyperplane with the unit normal v pointing into 
the positive halfspace. 

• Points. One considers points x G as hyperspheres of a vanishing 
radius. 

• Infinity. One compactifies the space by adding the point at infinity 
oo, with the understanding that a basis of open neighborhoods of 
oo is given, e.g., by the outer parts of the hyperspheres \x\'^ = r^. 
Topologically the so defined compactification is equivalent to a sphere 

• Contact elements. A contact element of a hypersurface is a pair con- 
sisting of a point x G M.^ and an (oriented) hyperplane P through x; 
alternatively, one can use a normal vector u to P at x. In the frame- 
work of Lie geometry, a contact element can be identified with a set (a 
pencil) of all hyperspheres S through x which are in an oriented con- 
tact with P (and with one another), thus sharing the normal vector v 
at X, see Fig. IHl 

All these elements are modelled in Lie geometry as points, resp. lines, 
in the (A^ -|- 2)-dimensional projective space P(R^+-^'^) with the space of 
homogeneous coordinates M^"*"^'^. The latter is the space spanned by A^ -|- 3 
linearly independent vectors ei,...,ejv+3 and equipped with the pseudo- 
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Figure 9: Contact element 

euclidean scalar product 

' 1, i = je{l,...,N + l}, 
{ei,ej) = l -1, i = je{N + 2,N + 3}, 

It is convenient to introduce two isotropic vectors 

^0 = — ^N+l), Goo = 5(eAr+2 + eAT+l), (3) 

for which 

(eo,eo) = (eoo,eoo) = 0, (eo,etx)) = 

The models of the above elements in the space M^+^'^ of homogeneous 
coordinates are as follows: 

• Oriented hypersphere with center c G and signed radius r E:M.: 

s = c + Bq + {\cf - r^)eoo + reN+s- (4) 

• Oriented hyperplane (v,x) = d with v € S^~^ and d G M: 

p = u + ■ eo + 2deoo +eAr+3- (5) 

• Point xeM.^: 

i; = a; + eo + [xpeoo + • eAr+3. (6) 
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• Infinity oo; 



do = Coo- 



(7) 



• Contact element {x, P): 

span{x,p). (8) 

In the projective space P(]R^+^'^) the first four types of elements are rep- 
resented by the points which are equivalence classes of (jlj)-© with respect 
to the relation ^ ~ <^ ^ = Ar/ with A € M* for ^,r/ G A con- 

tact element is represented by the line in P(R^+-^'^) through the points with 
the representatives x and p. We mention several fundamentally important 
features of this model: 

(i) All the above elements belong to the Lie quadric P(L^+"'^'^), where 

]L^+^'^ = {^eR''+''^:{tO=0}- (9) 

Moreover, points of P(L^+^'^) are in a one-to-one correspondence 
with oriented hyperspheres in R^, including degenerate case: proper 
hyperspheres correspond to points of P(L^"''^'-^) with both eg- and 
e^v+s-components non- vanishing, hyperplanes correspond to points of 
P(L''^'^^'^) with vanishing eo-component, points correspond to points 
of P(L^"*'"'^'^) with vanishing e7v+3-component, and infinity corresponds 
to the only point of P(L^"^^''^) with both eo- and ejv+3-components 
vanishing. 

(ii) Two oriented hyperspheres 81,82 are in an oriented contact (i.e., are 
tangent to each other with the unit normals at tangency pointing in 
the same direction), if and only if 

|ci -C2p = (n -r2)2, (10) 

and this is equivalent to (si, §2) = 0. 

(iii) An oriented hypersphere = {x G : |x — = r^} is in an oriented 
contact with an oriented hyperplane P = {x € : {v, x) = d}, if and 
only if 

{c,v)-r-d = Q. (11) 

Indeed, equation of the hyperplane P tangent to 5 at xq G 5* reads: 
(xo — c, X — c) = r^. Denoting hy v = {c — xo)/r the unit normal 
vector of P (recall that the positive radii are assigned to spheres with 
inward unit normals), we can write the above equation as {v,x) = d 
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with d = (c, (c — xo)/r) — r = {c,v) — r, which proves eq. ((TT|l . Now, 
the latter equation is equivalent to {s,p) = 0. 

(iv) A point x can be considered as a hypersphere of radius r = (in 
this case both oriented hyperspheres coincide). An incidence relation 
X € S with a hypersphere S (resp. x £ P with a hyperplane P) can 
be interpreted as a particular case of an oriented contact of a sphere 
of radius r = with S (resp. with P), and it takes place if and only if 
{x,s) = (resp. {x,p) = 0). 

(v) For any hyperplane P, there holds {6b, p) = 0. One can interpret 
hyperplanes as hyperspheres (of an infinite radius) through oo. More 
precisely, a hyperplane {v, x) = d can be interpreted as a limit, as 
r — >■ oo, of the hyperspheres of radii r with the centers located at 
c = rv + u, with {v, u) = d. Indeed, the representatives Q of such 
spheres are 

s = {rv + u) + + {2dr + {u, u))eoo + re^+s 

iv + 0(l/r)) + (l/r)eo + (2d + 0{l/r))e^ + e^+3 
= p + 0{l/r). 

Moreover, for similar reasons, the infinity 00 can be considered as a 
limiting position of any sequence of points x with |rE| — > 00. 

(vi) Any two hyperspheres ^i, ^2 in an oriented contact determine a con- 
tact element (their point of contact and their common tangent hyper- 
plane). For their representatives si, S2 in M^"*"^'^, the line in P(M^+^'^) 
through the corresponding points in P(L^''"^'^) is isotropic, i.e., lies en- 
tirely on the Lie quadric P(L^'^^'^). This follows from 

(aisi -I- 02^2, ai^i -I- 02^2) = 2aia2(si, S2) = 0. 

Such a line contains exactly one point whose representative x has van- 
ishing ejv+3-component (and corresponds to x, the common point of 
contact of all hyperspheres), and, if x 7^ 00, exactly one point whose 
representative p has vanishing eo-component (and corresponds to P, 
the common tangent hyperplane of all hyperspheres). In case when an 
isotropic line contains 60, all its points represent parallel hyperplanes, 
which constitute a contact element through 00. 

Thus, if one considers hyperplanes as hyperspheres of infinite radii, and 
points as hyperspheres of vanishing radii, then one can conclude that: 
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► oriented hyperspheres are in a one-to-one correspondence with points 
of the Lie quadric P(L^+-'^'^) in the projective space P(]R^+^'^); 

► oriented contact of two oriented hyperspheres corresponds to orthogo- 
naUty of (any) representatives of the corresponding points in P(R^+^'^). 

► contact elements of hypersurfaces are in a one-to-one correspondence 
with isotropic hnes in P(M'^'^^'^). We wiU denote the set of all such 
lines by . 

According to F. Klein's Erlangen Program, Lie geometry is the study of 
properties of transformations which map oriented hyperspheres (including 
points and hyperplanes) to oriented hyperspheres and, moreover, preserve 
the oriented contact of hypersphere pairs. In the projective model described 
above. Lie geometry is the study of projective transformations of P(]R^"''^'^) 
which leave P(L^^"^'^) invariant, and, moreover, preserve orthogonality of 
points of P(L^+^'^) (which is understood as orthogonality of their lifts to 
^ j^Af+1,2. gigg^j-iy^ relation does not depend on the choice of 
lifts). Such transformations are called Lie sphere transformations. 

Theorem 16 (Fundamental theorem of Lie geometry) 

a) The group of Lie sphere transformations is isomorphic to 0{N+1, 2)/{ib/}. 

b) Every line preserving dijjeomorphism of P(L^^^'^) is the restriction to 
P(L^"'"^'^) of a Lie sphere transformation. 

Since (non-) vanishing of the eg- or of the e7v+3-component of a point in 
P(L^"'"^'^) is not invariant under a general Lie sphere transformation, there 
is no distinction between oriented hyperspheres, oriented hyperplanes and 
points in Lie geometry. 

3.2 Mobius geometry 

Blaschke's book |Blj serves also as a classical source on Mobius geometry, a 
modern account can be found in jH Jj . 

Mobius geometry is a subgeometry of Lie geometry, with points distin- 
guishable among all hyperspheres as those of radius zero. Thus, Mobius 
geometry studies properties of hyperspheres invariant under the subgroup 
of Lie sphere transformations preserving the set of points. In the projec- 
tive model, points of are distinguished as points of P(L^^^'^) with the 
vanishing eAr+3-component. (Of course, one could replace here ejv+3 by 
any time-like vector.) Thus, Mobius geometry studies the subgroup of Lie 
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sphere transformations preserving the subset of P(L^^"'^'^) with the vanish- 
ing e/vT+s-component. Fohowing geometric objects in are elements of 
Mobius geometry. 

• (Non-oriented) hyperspheres S = {x G : |x — cp = r^} with centers 
c G and radii r > 0. 

• (Non-oriented) hyperplanes P = {x G M'^ : {v,x) = d}, with unit 
normals v € and d G M. 

• Points X G M^. 

• Infinity oo which compactifies into . 

In modelling these elements, one can use the Lie-geometric description 
and just omit the e^v+s-component. The resulting objects are points of the 
(A^ + l)-dimensional projective space P(M^"'"^'^) with the space of homoge- 
neous coordinates M^+-'^'^. The latter is the space spanned by + 2 linearly 
independent vectors ei, . . . ,ejv+2 and equipped with the Minkowski scalar 
product 

r 1, i = je{i,...,N+i}, 

{ei,ej) = I -1, i = j = N + 2, 
[ 0, i^j. 

We continue to use notations © in the context of the Mobius geometry. 
The above elements are modelled in the space M'^+^'^ of homogeneous co- 
ordinates as follows: 

• Hyper sphere with center c G and radius r > d: 

s = c + eo + (|cp-r2)eoo. (12) 

• Hyperplane {v, x) = d with v G S^^^ and d G M: 

p = w + • eo + 2deoo. (13) 

• Point X eM.^ : 

X = X + eo + jxpeoo- (14) 

• Infinity oo: 

6b = eoo. (15) 



26 



In the projective space P(M +^'-'^) these elements are represented by points 
which are equivalence classes of (fT^ - (fT^ with respect to the usual relation 
^ ~ <J4' ^ = Ar/ with A G M* for ^,r/ G M^+^'-^. Fundamental features of 
these identifications: 

(i) The infinity do can be considered as a limit of any sequence of x for 
X G with |x| — > oo. Elements x € U {00} are in a one-to-one 
correspondence with points of the projectivized light cone P(L^+^'^), 
where 

= {^gR^+i.i : (^,^) =0}. (16) 

Points X G correspond to points of P(L^^"^'"'^) with a non-vanishing 
eo-component, while 00 corresponds to the only point of P(L^+-'^'^) 
with the vanishing eQ-component. 

(ii) Hyperspheres s and hyperplanes p belong to P(IR^^"'^'^), where 

Cuf'' = UG]R^"'''':(e,O>0} (17) 

is the set of space-like vectors of the Minkowski space ]R^+^'-'^. Hy- 
perplanes can be interpreted as hyperspheres (of an infinite radius) 
through 00. 

(iii) Two hyperspheres Si , 52 with centers ci , C2 and radii ri , r2 intersect 
orthogonally, if and only if 

\ci - C2\^ = rf + rl (18) 

which is equivalent to (si,S2) = 0. Similarly, a hypersphere S inter- 
sects orthogonally with a hyperplane P, if and only if its center lies in 
P: 

{c,v)-d = 0, (19) 
which is equivalent to {s,p) =0. 

(iv) A point x can be considered as a limiting case of a hypersphere with 
radius r = 0. An incidence relation x (z S with a hypersphere S (resp. 
X G P with a hyperplane P) can be interpreted as a particular case 
of an orthogonal intersection of a sphere of radius r = with S (resp. 
with P), and it takes place if and only if (x, s) = (resp. {x,p) = 0). 

Note that a hypersphere S can also be interpreted as the set of points 
X G S*. Correspondingly, it admits, along with the representation s, the dual 
representation as a transversal intersection of P(L^^"^'^) with the projective 
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A?^-space P(s-'-), polar to the point s with respect to P(L^+^'-'^); here, of 
course, S"*- = {x G M^+^'^ : {s,x) = 0}. This can be generahzed to model 
lower-dimensional spheres. 

• Spheres. A /c-sphere is a (generic) intersection N — k hyperspheres Si 
{i = I, . . . , N — k). The intersection oi N — k hyperspheres represented 
by Si G M^^^'^ = 1, . . . , — /c) is generic if the {N — fc)-dimensional 
linear subspace of M^+^'-^ spanned by Sj is space-like: 

S = span(si, . . . , SM-k) C 

As a set of points, this A;-sphere is represented as P(L^+^'-'^nS-'-, where 

N~k 

= fl si = [xe M^+i'^ : (si, x) = . . . = {sN-k. x) = o} 

i=l 

is a (A;-|-2)-dimensional linear subspace of R^+^'^ of signature {k+1, 1). 

Through any k + 2 points xi, . . . , Xk+2 £ in general position one 
can draw a unique fe-sphere. It corresponds to the {k + 2)-dimensional 
linear subspace 

S-^ = span(fi, . . . ,Xjk+2), 

of signature (fc -|- 1, 1), with k + 2 linearly independent isotropic vec- 
tors xi,...,xa;+2 £ L^^^'^. In the polar formulation, this fc-sphere 
corresponds to the {N — A;)-dimensional space-like linear subspace 

k+2 

T,= [^xi = [se M^+^'^ : {s,xi) = ... = {s,Xk+2) = o}. 
1=1 

Mobius geometry is the study of properties of (non)-oriented hyper- 
spheres invariant with respect to projective transformations of P(R^"'"^'^) 
which map points to points, i.e., which leave P(L-^+^'^) invariant. Such 
transformations are called Mobius transformations. 

Theorem 17 (Fundamental theorem of Mobius geometry) 

a) The group of Mobius transformations is isomorphic to 0{N+1, 2± 
0~^{N + 1, 1), the group of Lorentz transformations of R^+^'^ preserving the 
time-like direction. 

b) Every conformal diffeomorphism of ~ U {oo} is induced by the 
restriction to P(L^'*"-'^'^) of a Mobius transformation. 
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The group 0^{N + 1,1) is generated by reflections 

As : M^+i'i ^ M^+i'i, Mx) = x- s. (20) 

{s,s) 

If s is a hypersphere (|12|1 . then the transformation induced on by 
is obtained from (^0]) by a computation with the representatives (fn)) for 
points and is given by: 

X I— > C + -j pr (x — c) (21) 

|x — 

(inversion in the hypersphere S = {x £ : |x — cp = r^}); similarly, if 
s = p is the hyperplane (|1,'-{|1 , then the transformation induced on by 
is easily computed to be 

2((v,x)-d) , , 

x^x- ^\ \ -V 22 



(reflection in the hyperplane P = {x G M : (t',x) = d}). 

Since (non-)vanishing of the eoo-component of a point in P(]R^"^^'-^) is 
not invariant under a general Mobius transformation, there is no distinction 
in Mobius geometry between hyperspheres and hyperplanes. 



3.3 Laguerre geometry 

Also in the case of Laguerre geometry the Blaschke's book |Blj serves as 
the indispensable classical source. One can find a modern account, e.g., in 

Laguerre geometry is a subgeometry of Lie geometry, with hyperplanes 
distinguishable among all hyperspheres, as hyperspheres through oo. Thus, 
Laguerre geometry studies properties of hyperspheres invariant under the 
subgroup of Lie sphere transformations which preserve the set of hyper- 
planes. Following objects in are elements of the Laguerre geometry. 

• (Oriented) hyperspheres S = {x £ : jx — = r^} with centers 
c G and signed radii r G M, can be put into correspondence with 
(iV + l)-tuples (c,r). 

• Points X G are considered as hyperspheres of radius zero, and are 
put into correspondence with (N + l)-tuples (x,0). 
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• (Oriented) hyperplanes P = {x G : {v,x) = d}, with unit normals 
V G and d E M., can be put into correspondence with {N + 1)- 
tuples {v,d). 

In the projective model of Lie geometry, hyperplanes are distinguished 

as elements of P(L^+-'^'^) with the vanishing eg-component. (Of course, one 
could replace here Bq by any isotropic vector.) Thus, Laguerre geometry 
studies the subgroup of Lie sphere transformations preserving the subset of 
P(L^+^'^) with the vanishing eo-component. 

There seems to exist no model of Laguerre geometry, where hypcrsphcres 
and hyperplanes would be modelled as points of one and the same space. De- 
pending on which of both types of elements is modelled by points, one comes 
to the Blaschke cylinder model or to the cyclographic model of Laguerre ge- 
ometry. We will use the first of these models which has an advantage of a 
simpler description of the distinguished objects of the Laguerre geometry, 
which are hyperplanes. The main advantage of the second model is a simpler 
description of the group of Laguerre transformations. 

The scene of the both models consists of two (iV-|- l)-dimensional projec- 
tive spaces, whose spaces of homogeneous coordinates, M^'^'^ and (M^'^'^)*, 
are dual to one another and arise from M^+^'^ by "forgetting" the Cq-, resp. 
Coo-components. Thus, M^'-*^'^ is spanned by AT -|- 2 linearly independent 
vectors ei,...,eiv, ejv-i-S) ^oo, and is equipped with a degenerate bilinear 
form of the signature (N, 1,1) in which the above vectors are pairwise or- 
thogonal, the first N being space-like: (ej,ej) = 1 for 1 < i < N , while 
the last two being time-like and isotropic, respectively: (gat+s, gat+s) = — 1 
and (eoo,eoo) = 0. Similarly, (M^'^'-^)* is assumed to have an orthogonal 
basis consisting of ei, . . . ,e7v, ^n+s, Gq, again with an isotropic last vector: 
(eo,eo) = 0. Note that one and the same symbol (•, •) is used to denote two 
degenerate bilinear forms in our two spaces. We will overload this symbol 
even more and use it also for the (non-degenerate) pairing between these 
two spaces, which is established by setting (eo,eoo) = additionally to 
the above relations. (Note that a degenerate bilinear form cannot be used 
to identify a vector space with its dual.) 

In both models mentioned above there holds: 

• Hyperplane P = {v, d) is modelled as a point in the space F{M.^'^'^) 
with a representative 

p = v + 2deoo + eiv+3- (23) 
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• Hypersphere S = (c, r) is modelled as a point in the space 
with a representative 

s = c + eo + reN+3. (24) 

Each one of the models appears if one considers one of the spaces as a pre- 
ferred (fundamental) one, and interprets the points of the second space as 
hyperplanes in the preferred one. In the Blaschke cylinder model, the pre- 
ferred space is the space P(]R''^'^'^) whose points model hyperplanes P C . 
A hypersphere S C is then modelled as a hyperplane E P(R^'^'-^) : 
(s,(^) = 0} in the space P(]R^'^'-^). Basic features of this model: 

(i) Oriented hyperplanes P C are in a one-to-one correspondence with 
points p of the quadric P(L^'^'-^), where 

l^N,i,i ^ 1^ ^ ^N,i,i . ^) = 0} . (25) 

(iii) Two oriented hyperplanes Pi,P2 C are in an oriented contact 
(parallel), if and only if their representatives pi, p2 differ by a vector 
parallel to eoo- 

(iii) An oriented hypersphere S C is in an oriented contact with an 
oriented hyperplane P C M^, if and only if p € s, that is, if {p, s) = 
0. Thus, a hypersphere S is interpreted as a set of all its tangent 
hyperplanes. 

The quadric P(L^'-^'^) is diffeomorphic to the Blaschke cylinder 

Z = {{v,d) G M^+^ : \v\ = 1} = S^"^ X M C R^+\ (26) 

Two points of this cylinder represent parallel hyperplanes, if they lie on one 
straight line generator of Z parallel to its axis. In the ambient space M^"*"^ 
of the Blaschke cylinder, oriented hyperspheres S C are in a one-to-one 
correspondence with hyperplanes non-parallel to the axis of Z: 

S ^ {{v,d) : {c,v) -d-r = 0}. (27) 

An intersection of such a hyperplane with 2. consists of points in Z which 
represent tangent hyperplanes to S" C M^, as follows from eq. 

In this paper, we will not use the cyclographic model of Laguerre geom- 
etry; its short description is put in Appendix El 
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4 Discrete curvature line parametrization in Lie, 
Mobius and Laguerre geometries 

Starting from here, we restrict ourselves to the geometry of surfaces in the 
three-dimensional Euclidean space M^. Accordingly, one should set iV = 3 
in all previous considerations. 

It is natural to consider following objects as discrete surfaces in the 
various geometries discussed above: 

• In Lie geometry, a surface is viewed as built of its contact elements. 
These contact elements are interpreted as points of the surface and 
tangent planes (or, equivalently, normals) at these points. This can be 
discretized in a natural way: a discrete surface is a map 

{x,P) : 1? ^ {contact elements of surfaces in M^}, 

or, in the projective model of Lie geometry, a map 

^ (28) 

where, recall, £q'^ denotes the set of isotropic lines in P(M^'^). 

• In Mobius geometry, a surface is viewed simply as built of points. A 
discrete surface is a map 

x: ^ R^, 
or, in the projective model, a map 

X : ^ P(L4'1). (29) 

• In Laguerre geometry, a surface is viewed as the envelope of the system 
of its tangent planes. A discrete surface is a map 

P : 1? ^ {oriented planes in R^}, 

or, in the projective model, a map 

^ P(L3.i.1). (30) 

It should be mentioned that a substantial part of the description of a 
surface in Laguerre geometry is its Gauss map 

v.I?^ S^, (31) 

consisting of unit normals v to the tangent planes P = {v,d). 

Thus, description of a discrete surface in Lie geometry contains more 
information than description of a discrete surface in Mobius or in Laguerre 
geometry. Actually, the former merges the two latter ones. 
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4.1 Lie geometry 



The following definition is a discretization of the Lie-geometric description 
of curvature line parametrized surfaces, as found, e.g., in |Blj . 

Definition 18 (Principal contact element nets. Euclidean model) 

A map 

(x, P) : — > {contact elements of surfaces in M^} 

is called a principal contact element net, if any two neighboring contact 
elements {x,P), {xi,Pi) have a sphere S^^^ in common, that is, a sphere 
touching both planes P, Pi at the corresponding points x, Xi. 

Thus, the normals to the neighboring planes P, Pi at the corresponding 
points X, Xi intersect at a point c^*-* (the center of the sphere S"*-*^), and the 
distances from c^*) to x and to Xi are equal, see Fig. The spheres S*^*) , 
attached to the edges of 1? parallel to the z-th coordinate axis, will be called 
principal curvature spheres of the discrete surface. 




Figure 10: Principal curvature sphere 

A direct translation of Definition into the projective model looks as 
follows: 

Definition 19 (Principal contact element nets. Projective model) 

A map i : W? ^ ^Lq'^ is called a principal contact element net, if it is a 
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discrete congruence of isotropic lines inP(IR^'^), that is, any two neighboring 
lines intersect: 

i{u)n£{u + ei) = s^'Hu) e¥{h'^'^), VuGZ^, Vi = l,2. (32) 

In the projective model, the representatives of the principal curvature 
spheres S*^*-* of the i-th coordinate direction build the corresponding focal 
net of the line congruence i; 

: ^2 ^ P(l4.2)^ i = l,2, (33) 

cf. Definition 1111 According to Theorem 1121 both focal nets are Q-nets in 
P(M^'^). This motivates the following definition. 

Definition 20 (Discrete R-congruence of spheres) A map 

S : {oriented spheres in M^} 

is called a discrete R-congruence (Rihaucour congruence) of spheres, if the 
corresponding map 

s : ^ P(L4.2) 

is a Q-net in P(M^'2). 

A geometric characterization of discrete R-congruences will be given in Sect. 

m 

Corollary 21 (Curvature spheres build an R-congruence) For a dis- 
crete contact element net, the principal curvature spheres of the i-th coordi- 
nate direction (i=l,2) build a two-dimensional discrete R-congruence. 

Turning to transformations of principal contact element nets, we intro- 
duce the following definition. 

Definition 22 (Ribaucour transformation. Euclidean model) Two 

principal contact element nets 

{x,P), (x~^,P~^) : — > {contact elements of surfaces in M'^} 

are called Ribaucour transforms of one another, if any two corresponding 
contact elements {x,P) and {x~^, P^) have a sphere S in common, that is, a 
sphere which touches both planes P, P~^ at the corresponding points x, x~^ . 

Again, a direct translation of Definition EH into the language of the 
projective model gives: 
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Figure 11: Ribaucour transformation 



Definition 23 (Ribaucour transformation. Projective model) Two 

principal contact element nets £, : 1? ^ XLq^ are called Ribaucour trans- 
forms of one another, if these discrete congruences of isotropic lines are in 
the relation of F-transformation, that is, if any pair of the corresponding 
lines intersect: 

l{u) n = s(n) e P(L^'2), Vu € I?. (34) 

Spheres of a Ribaucour transformation are attached to the vertices u 
of the lattice 1?, or, better, to the "vertical" edges connecting the vertices 
(u, 0) and (ti, 1) of the lattice J? x {0,1}. In the projective model, their 
representatives 

s:I? ^ P(L'^'2) (35) 

build the focal net of the three-dimensional line congruence for the third 
coordinate direction. From Theorem 1121 there follows: 

Corollary 24 (Spheres of a Ribaucour transformation build an R- 
congruence) The spheres of a generic Ribaucour transformation build a 
discrete R-congruence. 

Now, we turn to the study of the geometry of an elementary quadrilateral 
of contact elements of a principal contact element net, consisting of £ ~ 
(x,P), £i ~ (xi,Pi), £2 ~ {X2,P2), and £12 ~ (xi2,Pi2). 
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We leave aside a degenerate umbilic situation, when all four lines have 
a common point and span a four-dimensional space. Geometrically, this 
means that one is dealing with four contact elements of a sphere S* C M^. In 
this situation, one cannot draw any further conclusion about the four points 
j;, j;2) ^12 on the sphere S: they can be arbitrary. 

In the non-umbilic situation, the space spanned by the four lines i, ^i, ^2, ^12 
is three-dimensional. The four elements ^i, ^^12 £ P(L^'^) correspond- 
ing to the points x, xi, X2, 2:12 S are obtained as the intersection of 
the four isotropic lines ^,^1,^21^12 with the projective hyperplane P(e^) 
in P(M^'-^). Therefore, the four elements x,xi,X2,a:i2 lie in a plane. A suit- 
able framework for the study of this configuration is the projective model 
of the Mobius geometry. Namely, omitting the inessential (vanishing) eg- 
component, we arrive at a planar quadrilateral in the Mobius sphere P(L^'^). 
We devote Sect. 14.21 to the study of such objects. 

Analogously, the four elements p,Pi,P2,Pi2 £ P(L'^'^) corresponding to 
the planes P,Pi,P2,Pi2 G P^ are obtained as the intersection of the four 
isotropic lines ^,^1,^2,^12 with the projective hyperplane P(e;^) in P(R^'^). 
Therefore, also the four elements p,Pi,P2,Pi2 lie in a plane. A suitable 
framework for the study of such a configuration is the projective model of 
the Laguerre geometry; this will be performed in Sect. 14.31 

4.2 Mobius geometry: circular nets 

Circular nets were introduced and studied in the context of integrable sys- 
tems in p iCDSllKSj . 

Caution: in this section, the notation x refers to the Mobius-geometric 
representatives in L^'^, and not to the Lie-geometric ones in L'''^. The former 
are obtained from the latter one by omitting the (vanishing) eg-component. 

We assume that the principal contact element nets under consideration 
are generic, i.e., do not contain umbilic quadruples. The main result of this 
section is the following claim. 

Theorem 25 (Points of principal contact element nets form circu- 
lar nets) For a principal contact element net 

{x,P) : 1? {contact elements of surfaces in M^}, 

its points X : ^ P^ form a circular net. 

This statement refers to the notion, which can be defined in two different 
ways. 
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Definition 26 (Circular net. Euclidean model) A net x : is 

called circular, if the vertices of any elementary quadrilateral {x,Xi,Xij,Xj) 
(at any u € Z"* and for all pairs 1 < i ^ j < m) lie on a circle (in particular, 
are co-planar). 

Definition 27 (Circular net. Projective model) A net x : Z*" — > is 

called circular, if the corresponding x : Z™ — > P(L ' ) is a Q-net in P(M^'i). 

This time a translation between the Euclidean model and the projective 
model is not straightforward and actually constitutes the matter of Theorem 
1251 indeed, this theorem has already been demonstrated (or, better, is ob- 
vious) in terms of Definition 1271 and it remains to establish the equivalence 
of Definitions EHEZl 

Conceptual proof. The linear subspace of M^'^ spanned by the isotropic 
vectors three-dimensional. Its orthogonal complement is 

therefore two-dimensional and lies in M^^^. Therefore, it represents a circle 
(an intersection of two spheres). □ 

Computational proof. For arbitrary representatives x G L^'^ of x, the 
requirement of Definition |2Z1 is equivalent to equation of the type (CQ) . Since 
the representatives x = x-t-eo+jxpCoo fixed in ()14() lie in an affine hyperplane 
of M'*'^ (their eo-component is equal to 1), one has an equation of the type 
((21) for them. Clearly, this holds if and only if x is a Q-net in and |xp 
satisfies the same equation ^ as x does. We show that the latter condition 
is equivalent to circularity. On a single planar elementary quadrilateral 
{x,Xi,Xij,Xj), the function satisfies eq. ((2j simultaneously with \x — 
cp = |xp — 2{x, c) + |cp with any c G R^. Choose c to be the center of the 
circle through the three points X, Xj, Xj, so that |x — = |xj — cp = |xj — cp. 
Then eq. ^ for |x — cp turns into \xij — cp = |x — cp, which means that 
Xij lies on the same circle. □ 

Two-dimensional circular nets (m = 2) are discrete analogs of the curva- 
ture lines parametrized surfaces, while the case m = 3 discretizes orthogonal 
coordinate systems in R^. A construction of an elementary hexahedron of a 
circular net is based on the following geometric theorem: 

Theorem 28 (Elementary hexahedron of a circular net) Given seven 
points X, Xi, and Xij (1 < « < i < 3) in R'^, such that each of the three 
quadruples {x,Xi,Xj,Xij) lies on a circle Cij, define three new circles TiCjk 
as those passing through the triples {xi,Xij,Xik), respectively. Then these 
new circles intersect at one point, see Fig. fTH - 

2^123 = nC23 n T2C31 n T3C12 . 
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Figure 12: An elementary hexahedron of a circular net 



Proof. This is a particular case of Theorem El applied to the quadric 
P(L4'1). □ 

This theorem can be proven also by elementary geometrical considera- 
tions. If one notes that under conditions of Theorem EH] the seven points 
X, Xi, Xij lie on a two-dimensional sphere, and performs a stereographic 
projection of this sphere with the pole at x, one arrives at a planar picture 
which is nothing but the classical Miquel theorem. 

4.3 Laguerre geometry: conical nets 

Conical meshes have been introduced recently in |LPWYW] . 

Caution: in this section, the notation p refers to the Laguerre-geometric 
representatives in L^'^'^, and not to the Lie-geometric ones in L"^'^. The for- 
mer are obtained from the latter by omitting the (vanishing) eo-component. 

As in the previous section, we assume that the principal contact element 
nets under consideration do not contain umbilic quadruples. The main result 
of this section is the following claim. 

Theorem 29 (Tangent planes of principal contact element nets 
form conical nets) For a principal contact element net 

{x,P) : Tl? {contact elements of surfaces in M^}, 

its tangent planes P : Z™" —>■ {oriented planes in M^} form a conical net. 

This statement refers to the notion which can be defined in two different 
ways. 
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Definition 30 (Conical net. Euclidean model) A net P : — > 

{oriented planes in M"^} is called conical, if at any u S and for all 
pairs 1 < i j < m the four planes P, Pi, Pij, Pj touch a cone of revolution 
(in particular, intersect at the tip of the cone). 

Definition 31 (Conical net. Projective model) A net P : 17^ — s- 

{oriented planes in M'^} is called conical, if the corresponding p : Z™ —>■ 
P(L^'i'^) is a Q-net in P(M3'1'1). 

Theorem l29l is obvious in terms of Definition l311 so the real content of this 
theorem is the translation between the Euclidean model and the projective 
model, that is, establishing equivalence of Definitions l3(H I^Tl 

Proof. Representatives p in (|23|) form a Q-net, if and only if they 
satisfy eq. ©, that is, if v : and d : Z'" M satisfy this equation. 

Equation Q for v yields that v : Z™" — > §^ is actually a Q-net in S^, so 
that any quadrilateral {v,Vi,Vij,Vj) in is planar and therefore circular. 
Equation @ for {v,d) yields that the (unique) intersection point of three 
planes P, Pi, Pj lies on Pij, as well, so that all four planes intersect in 
one point. Thus, we arrived at a characterization of conical nets in the 
sense of Definition |^ as those nets of planes for which every quadruple of 
planes {P, Pi, Pij, Pj) is concurrent and every quadrilateral {v,Vi,Vij,Vj ) of 
unit normal vectors is planar. It is clear that this description is equivalent to 
that of Definition 1301 The direction of the axis of the tangent cone coincides 
with the spherical center of the quadrilateral {v,Vi,Vij,Vj) in 

Thus, conical nets are Q-nets with circular Gauss maps. It is worthwhile 
to mention that, in order to prescribe a conical net, it is enough to prescribe 
a circular Gauss map v : Z™ — > and additionally the numbers d (i.e., the 
planes P = {v,d)) along the coordinate axes of Z™. Indeed, these data allow 
one to reconstruct the conical net uniquely. This is done via a recursive 
procedure, whose elementary step consists in finding the fourth plane Pij 
provided three planes P,Pi,Pj and the normal direction Vij of the fourth 
one are known. But this is easy: Pij is the plane normal to Vij through the 
unique intersection point of the three planes P, Pi , Pj . 

4.4 Synthesis 

In view of Theorems 1251 1291 it is natural to ask whether, given a circular 
net X : Z^ ^ M^ or a conical net P : Z^ ^ {oriented planes in M^}, there 
exists a principal contact element net 

{x,P) : 1? {contact elements of surfaces in M^}, 
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with the prescribed half of the data (x or P). A positive answer to this 
question is a corollary of the following general theorem. 

Theorem 32 (Extending R-congruences of spheres to curvature 
line parametrized surfaces) Given a discrete R-congruence of spheres 

S :W? ^ {oriented spheres in M'^}, 

there exists a two-parameter family of principal contact element nets 

(x,P) : 1? — > {contact elements of surfaces in M^} 

such that S belongs to the contact element {x,P), i.e., P is the tangent 
plane to S at the point x E S", for all n G Z^. Such a principal contact ele- 
ment net is uniquely determined by prescribing a contact element {x, P){0, 0) 
containing the sphere S{0,0). 

Proof. The input data is a Q-net s : Z,"^ —>■ P(L^'^) in the Lie quadric, 
and we are looking for a congruence of isotropic lines t : J? ^ iLg'^ such 
that s{u) G i{u) for all u The construction starts with an arbitrary 

isotropic line ^(0,0) through s(0,0), and hinges on the following lemma. 

Lemma 33 For an isotropic line i G -Cq'^ and a point si G P(L'*'^) not lying 
on i, there is a unique isotropic line £i through si intersecting i. 

Proof. Let s,a be two arbitrary points on i (in homogeneous coordinates), 
so that the line i is given by the linear combinations as + Pa. Relation 

(as + f3a, si) = yields 

a : (3 = - {a,si) : {s, si). 

Thus, there exists a unique point s^^^ G i such that {s^^\ si) = 0. Now £i is 
the line through si and s^^\ □ 

Proof of Theorem I32L continued. With the help of Lemma ESI one can 
construct the isotropic lines of the congruence along the coordinate axes, 

£:Zx {0} ^ £^'2 and £ : {0} x Z /Lg'^. 

Next, one has to extend the congruence £ from the coordinate axes to the 
whole of 7?. An elementary step of this extension consists in finding, for 
three given isotropic lines £, £i, £2 (such that £ intersects both £1 and £2) 
the fourth one, £12, intersecting £1 and £2 and going through a given point 
512- One can use for this Lemma but then one has to demonstrate 
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that this construction is consistent, i.e., that the hnes ii2 obtained from 
the requirements of intersecting with ii and with £2 coincide. We show 
this with the fohowing argument. The space V = span(^, £1, £2) is three- 
dimensional. The points s, si, S2 he in V. By the hypothesis of the theorem, 
the quadrilateral (s, si, S12, S2) is planar, therefore §12 lies in V, as well. 
Draw two planes in V: Hi = span(^i,si2) and 112 = span(£2, •S12). Their 
intersection is a line £12 through S12. It remains to prove that this line is 
isotropic. For this, note that ii2 can be alternatively described as the line 
through two points s^^^ = £1 H lyi and s~2^ = ^2 H l\2- Both these points lie 
in P(L^'^), since they belong to the isotropic lines l\ and £2^ respectively. 
But it is easy to see that a line in P(R^'^) through two points from P(L^'-^) 
is either isotropic, or contains no further points from P(L^'^), depending on 
whether these two points are polar to one another (with respect to P(L'^'^)) 
or not. In our case the line i\2 contains, by construction, one further point 
S12 from P(L^'^), therefore it has to be isotropic. □ 

Since the representatives x in P(L^'^) of a circular net x : ^ form 
a Q-net in P(R^'^), and the same holds for the representatives j5 in P(L^'-^) 
of a conical net P : 7? ^ {oriented planes in M^}, we come to the following 
conclusion (obtained independently by Pottmann ^). 

Corollary 34 (Extending circular and conical nets to principal con- 
tact element nets) 

i) Given a circular net x : — > , there exists a two-parameter family 
of conical nets P : — > {planes in M^} such that x G P for all u £ 7?, and 
the contact element net 

{x,P) : 1? — > {contact elements of surfaces in M^} 

is 'principal. Such a conical net is uniquely determined by prescribing a plane 
P(0, 0) through the point x(0, 0). 

ii) Given a conical net P : 1? ^ {oriented planes in M^}, there exists a 
two-parameter family of circular nets x : — > such that x G P for all 
u S 1? , and the contact element net 

(x,P) : 1? — > {contact elements of surfaces in R^} 

is principal. Such a circular net is uniquely determined by prescribing a 
point x(0,0) in the plane P(0,0). 

These relations can be summarized as in Fig. 1131 Note that the axes of 
conical nets corresponding to a given circular net coincide with the Gauss 
map at its vertices, considered by Schief |Sch2j . 



41 




LAGUERRE 
LIE 

MOBIUS 



Figure 13: Elementary quadrilateral of a curvature line parametrized sur- 
face with vertices x and tangent planes P in the projective model. The 
vertices x build a circular net (Mobius geometry), and lie in the planes P 
building a conical net (Laguerre geometry). Contact elements (.x, P) are 
represented by isotropic lines I (Lie geometry). Principal curvature spheres 
S'(') pass through pairs of neighboring points x, Xi and are tangent to the 
corresponding pairs of planes P, Pi. 
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Remark. In the situations of Corollary 1341 i.e., when the R-congruence 
S consists of points x (and is therefore a circular net) or of planes P (and 
is therefore a conical net), the elementary construction step of Lemma 15^ 
allows for a very simple description from the Euclidean perspective in W^. 
This has been given by Pottmann 0. 

i) Given a contact element (x, P) and a point xi, find a plane Pi through 
xi so that there exists a sphere S*-"*^-* tangent to both planes P, Pi at 
the points x, xi, respectively. Solution: Pi is obtained from P by the 
reflection in the bisecting orthogonal plane of the edge The 
center c^^^ of the sphere S^^^ is found as the intersection of the normal 
to P at X with the bisecting orthogonal plane of the edge [x,xi]. 

ii) Given a contact element {x,P) and a plane Pi, find a point xi in Pi 
so that there exists a sphere S^^^ tangent to both planes P, Pi at the 
points X, xi, respectively. Solution: the point xi is obtained from x 
by the reflection in the bisecting plane of the dihedral angle formed by 
P, Pi. The center c^^^ of the sphere S^^^ is found as the intersection 
of the normal to P at x with this bisecting plane. 

5 R-congruences of spheres 

In Sect. |1J Corollaries 1211 OH we have seen that principal curvature spheres 
of a principal contact element net and spheres of a Ribaucour transformation 
build discrete R-congruences, introduced in Definition EOl In this section we 
study the geometry of discrete R-congruences of spheres. Definition [20l can 
be re-formulated as follows: a map 

S '.l? ^ {oriented spheres in M^}, 

or the corresponding map 

into the space of homogeneous coordinates, is called a discrete R-congruence 
of spheres, if for any u £ Z"* and for any pair 1 < i ^ j < m the linear 
subspace 

S = span{s, Si, Sj,Sij) 

is three-dimensional. Thus, to any elementary square of Z*^ there corre- 
sponds a three-dimensional linear subspace S C M^'^. 
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The R-congruence of principal curvature spheres 5^*^ of the i-th coordi- 
nate direction is degenerate in the sense that the subspaces of its elementary 
quadrilaterals 

S = span(s»,s«,s«,sf) 

contain two-dimensional isotropic subspaces (corresponding to ii and iij). 
The R-congruence of spheres of a generic Ribaucour transformation is, on the 
contrary, non-degenerate: its S's do not contain two-dimensional isotropic 
subspaces, and its elementary quadrilaterals are included in planar families 
of spheres, introduced in the following definition. 

Definition 35 (Planar family of spheres) A planar family of spheres 

is a set of spheres whose representatives s € P(L'^'^) are contained in a 
projective plane P(S), where Yj is a three-dimensional linear subspace of 
M^'^ such that the restriction of (•, •) to S is non- degenerate. 

Thus, a planar family of spheres is an intersection P(EnL^'^). Clearly, there 
are two possibilities: 

• Signature of (•, ■)]■£ is (2, 1), so that signature of (•, ■)\j^± is also (2, 1). 

• Signature of (•, ■)]■£ is (1, 2), so that signature of (•, ■)\-s± is (3, 0). 

It is easy to see that a planar family is one-parametric, parametrized by a 
circle S^. Indeed, if ei, 62, 63 is an orthogonal basis of E such that (ei, ei) = 
(62,62) = —(63,63) = 1 (say), then the spheres of the planar family come 
from the linear combinations s = aiei + 0262 + 63 with 

{aiei + 0262 -I- 63, CKiei -|- 0:262 -|- 63) = <^ a\-\- a\ = l. 

In the second of the cases mentioned above, the space S"*- only has a trivial 
intersection with L^'^, so that the spheres of the planar family P(L^'^ n E) 
have no common touching spheres. This case has no counterpart in the 
smooth differential geometry. From the point of view of discrete differential 
geometry the first case is more significant. 

Definition 36 (Cyclidic family of spheres) A planar family of spheres 
is called cyclidic, if the signature of (•, •)|s is (2, 1), so that the signature of 
(•,-)ls-L is also (2,1). 

Thus, for a cyclidic family P(L'^'^ n S) there is a dual cyclidic family 
P(L^'^ n S-*-) such that any sphere of the first one is in an oriented contact 
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with any sphere of the second one. The family P(L^'^ H S), as any one- 
parameter family of spheres, envelopes a canal surface in M^, and this surface 
is an envelope of the dual family P(L^'^ n S"*-), as well. Such surfaces are 
called Dupin cyclides. Thus, to any elementary quadrilateral of a discrete 
R-congruence whose spheres {s, Si, Sij, Sj) span a subspace of the signature 
(2,1) there corresponds a Dupin cyclide. 



Figure 14: A cyclidic family of spheres through a circle 
Examples: 

► points of a circle build a planar cyclidic family of spheres (of radius 
zero). The dual family consists of all (oriented) spheres through this circle, 
with centers lying on the line through the center of the circle orthogonal 
to its plane, see. Fig. I14L The corresponding Dupin cyclide is the circle 
itself. It can be shown that any Dupin cyclide is an image of this case under 
a Lie sphere transformation. For a circular net, considered as a discrete 
R-congruence, each elementary quadrilateral carries such a structure. 

► planes tangent to a cone of revolution build a planar cyclidic family 
of spheres, as well. The dual family consists of all (oriented) spheres tangent 
to the cone, with centers lying on the axis of the cone, see. Fig. 1151 The 
corresponding Dupin cyclide is the cone itself. For a conical net, considered 
as a discrete R-congruence, each elementary quadrilateral carries such a 
structure. 

Theorem 37 (Common tangent spheres of two neighboring quadri- 
laterals of an R-congruence) For two neighboring quadrilaterals of a 
discrete R-congruence of spheres, carrying cyclidic familes, there are gener- 
ically exactly two spheres tangent to all six spheres of the congruence. 

Proof. Let the quadrilaterals in question belong to the planar families gen- 
erated by the subspaces Si and T,2 of the signature (2,1). These quadrilat- 
erals share two spheres si and ,§2, which span a linear space of the signature 
(1,1). Each of the planar families Si and S2 adds one space-like vector, 
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Figure 15: A cyclidic family of spheres tangent to a cone 

so that the hnear space Si U S2 spanned by all six spheres of the congru- 
ence is four-dimensional and has the signature (3,1), so that its orthogonal 
complement (Si U S2)"'" is two-dimensional and has the signature (1,1). In- 
tersection of L^'^ with a two-dimensional linear subspace of the signature 
(1,1) gives, upon projectivization, exactly two spheres: indeed, if 61,62 form 
an orthogonal basis of (Si U S2)"'" with (61,61) = —(62,62) = 1, then the 
spheres in this space correspond to ai6i + 0262 with 

{aiei + 0262,0161 + 0262) = <^ of = 0.2 <^ ai : a2 = il- D 
In particular: 

► For any two neighboring quadrilaterals of a circular net, there is one 
non-oriented sphere (hence two oriented spheres) containing both circles. Its 
center is the intersection point of the lines passing through the centers of 
the circles orthogonallly to their respective planes, see Fig. 1 161 

► For any two neighboring quadrilaterals of a conical net, there is a 
unique oriented sphere touching both cones (the second such sphere is the 
point at infinity). The center of this sphere is the intersection point of the 
axes of the cones. 

The next theorem is proved in exactly the same way as Theorem 1371 

Theorem 38 (Common tangent spheres of an elementary hexa- 
hedron of an R-congruence) For an elementary hexahedron of a discrete 
R-congruence of spheres, with all faces carrying cyclidic families, there are 
generically exactly two spheres tangent to all eight spheres at its vertices. 

It should be mentioned that these spheres, attached to elementary hex- 
ahedra, do not form a discrete R-congruence, contrary to what has been 
alleged by A. Doliwa as a main result of |D3j . 
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Figure 16: Normals of two neighboring quadrilaterals of a circular net in- 
tersect: both lie in the bisecting orthogonal plane of the common edge 




Figure 17: Normals of two neighboring quadrilaterals of a conical net inter- 
sect: two common planes of the quadrilaterals arc tangent to both cones, 
therefore the axes of both cones lie in the bisecting plane of these two planes 
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We now turn to a geometric characterization of discrete R-congruences. 
From eq. ^ there fohows immediately that a map 

S : {oriented spheres in M^}, 

is a discrete R-congruence, if and only if the centers c : M'^ form a 

Q-net in R^, and the two real- valued functions, 

|cp -r^ :Z'^ and r : Z"" ^ M, 

satisfy the same equation of the type © as the centers c. By omitting 
the latter requirement for the signed radii r, one comes to a less restrictive 
definition than that of R-congruence. Actually, this definition belongs to 
Mobius geometry and refers to notations of Sect. 13.21 (with = 3). 

Definition 39 (Q-congruence of spheres) A map 

S : Z™" —)■ I non-oriented spheres in M^}, (36) 
is called an Q-congruence of spheres, if the corresponding map 

s:Z^^ P(R^;\), s = c + eo+ (|cp - r^)eoo, (37) 
is a Q-net in P(R^'^). 

Thus, a map (|36|) is a Q-congruence, if and only if the centers c : — > 
of the spheres S form a Q-net in R^, and the function |cp — satisfies the 
same equation Q as the centers c. 

Theorem 40 (Characterization of R- among Q-congruences) Four 
(oriented) spheres {S, Si, Sij, Sj) in R^ comprise an elementary quadrilateral 
of an R-congruence, if and only if they comprise (as non-oriented spheres) 
an elementary quadrilateral of a Q-congruence, and satisfy additionally the 
following condition: 

(R) There exists a non-point sphere in an oriented contact with all four 
oriented spheres S, Si, Sj, Sij. 

Under this condition, any sphere in an oriented contact with three spheres 
S,Si,Sj is in an oriented contact with the fourth one, Sij, as well. 
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Proof. Let Sq be a sphere with the center cq and (finite) oriented radius 
ro 7^ in an oriented contact with the three spheres S,Si,Sj. This means 
that the following conditions are satisfied: 



(c,co)-i(|cp-r2)-i(|co 



ro) - rro = 



(38) 



(tangency of S, Sq, cf. (|Tn|) l. and two similar equation with (c, r) replaced 
by {ci,ri) and {cj,rj). Now, using the fact that c and |cp — satisfy one 
and the same equation of the type ((21), we conclude that eq. (j^H]) is fulfilled 
for {cij,rij), if and only if r satisfies the same equation as c and |cp — 
do. This proves the theorem in the case when the common tangent sphere 
5o for the three spheres 5, Si, Sj has a finite radius. The case when Sq has 
an infinite radius, i.e., is actually a plane, is dealt with similarly, with the 
help of equation 



which comes to replace eq. (|HH|) . □ 

Remark. We have already seen that, generically, if three oriented 
spheres S, Si, Sj have a common sphere in an oriented contact, then they 
have a one-parameter (cyclidic) family of common touching spheres, repre- 
sented by a three-dimensional linear subspace S of R^'^. It is easy to see 
that if the projection of S onto is non- vanishing, then the family of 
spheres represented by S-*- contains exactly two planes. (The only excep- 
tional case is that of a conical cyclidic family S, all of whose elements have 
vanishing eo-component and represent planes, while the family contains 
no planes.) Therefore, in all cases but the conical one, condition (R) can be 
replaced by the following requirement: 

(Ro) The four oriented spheres S, Si, Sj, Sij have a common tangent plane 
(actually, two common tangent planes). 

It remains to give a geometric characterization of Q-congruences. This 
is done in the following theorem. 

Theorem 41 (Three types of Q-congruences) Four (non- oriented) 
spheres {S, Si, Sij, Sj) in comprise an elementary quadrilateral of a Q- 
congruence, if and only if they satisfy one of the following three conditions: 

(i) they have a common orthogonal circle, or 

(a) they intersect along a pair of points (a 0-sphere), or else 

(Hi) they intersect at exactly one point. 




(39) 
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Case (in) can be regarded as a degenerate case of both (i) and (ii). 

Conceptual proof. Caution: notations in this proof refer to the Mobius- 
geometric objects which are different from the Lie-geometric objects denoted 
by the same symbols. The hnear subspace S of M^'"*^ spanned by the points s, 
Si, Sj, Sij is three-dimensional, so that its orthogonal complement S-*- is two- 
dimensional. If S"*- lies in MqJj\, i.e., if the restriction of the Minkowski scalar 
product to S"*- is positive-definite (of signature (2,0)), then S"*- represents 
a 1-sphere (a circle) orthogonal to our four spheres, and we have the case 
(i). If, on the contrary, the restriction of the scalar product to has 
signature (1, 1), so that S lies in Mq^\, then S represents a 0-sphere which 
is the intersection of our four spheres, and we have the case (ii). Finally, 
if the restriction of the scalar product to S is degenerate, then S n S-*- is 
an isotropic one-dimensional linear subspace, which represents the common 
point of our four spheres, and we have the case (iii). □ 

Computational proof. The quadrilateral in with the vertices at the 
sphere centers c, Cj, cj, Cij is planar; denote its plane by II. In the same way 
as in the proof of Theorem 23 we show that there is a point C G II such 
that 

|c - C|2 -r^ = \ci - - rf = \cj - C\^ - r| = \cij - Cf - r%. (40) 

Indeed, the first two of these equations define C uniquely as the intersection 
of two lines ii and in II, where 

■^i = {x G n : (2x — Ci — c,Ci — c) = r"^ — r^^}, 

and then the last equation in (|40j) is automatically satisfied. If the common 
value of all four expressions in (|in|) is positive (say, equal to i?^), then the 
four spheres under consideration are orthogonal to the circle in the plane II 
with the center C and radius i?, so that we have the case (i), see Fig. 1181 
If the common value of ()40|) is negative (say, equal to — i?^), then the pair 
of points on the line through C orthogonal to II, at the distance R from 
C, belong to all four spheres, so that we have the case (ii). Finally, if the 
common value of (|40)) is equal to 0, then C is the intersection point of all 
four spheres, and we have the case (iii). □. 

Clearly, case (i) of Q-congruences reduces to circular nets, if the radii 
of all spheres become infinitely small, cf. Fig. 1181 Q-congruences with 
intersections of type (ii) are natural discrete analogs of sphere congruences 
parametrized along principal lines. 



50 



Figure 18: Elementary quadrilateral of a Q-congruence of spheres, the or- 
thogonal circle case 

Some remarks about Q-congruences of spheres are in order here. They 
are multidimensionally consistent, with the following reservation: given 
seven points s, Si, Sij in P(Mq^^), the Q-property (planarity condition) 
uniquely defines the eighth point S123 in P(M^'^), which, however, might get 
outside of P(RQy^-), and therefore might not represent a real sphere. Thus, 
the corresponding discrete 3D system is well-defined on an open subset of the 
space of initial data only. As long as it is defined, it can be used to produce 
transformations of Q-congruences, with usual pcrmutability properties. 

Note the following difference between Q-congruences and R-congruences: 
given three spheres S, Si, Sj of an elementary quadrilateral, one has a two- 
parameter family for the fourth sphere Sij in the case of a Q-congruence, 
and only a one-parameter family in the case of an R-congruence. This is 
a consequence of the fact that Mq^\ is an open set in M^'^, while L^'^ is a 
hypersurface in M^'^. 

A Appendix: cyclographic model of Laguerre 
geometry 

In the cyclographic model of Laguerre geometry, the preferred space is the 
space of hyperspheres (M^'^'^)*, so hyperspheres S C M-^ are modelled as 
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points s G P((M-^'^'-'^)*), while hyperplanes P C are modelled as hyper- 
planes : (j5,^) = 0} C P((R^'^'^)*). Thus, a hypcrplanc P is interpreted 
as a set of hypersphercs S which are in oriented contact with P. 
Basic features of this model: 

(i) The set of oriented hyperspheres S C is in a one-to-one correspon- 
dence with points 

a = {c,r) (41) 

of the Minkowski space M^'-^ spanned by the vectors ei, . . . , e^v, bn+s- 
This space has an interpretation of an afhne part of P((R^'^'^)*) . 

(ii) Oriented hyperplanes P C can be modelled as hyperplanes in R^'^: 

TT = {(c,r) € R^'^ : ((w, 1), (c, r)) = {v,c) - r = d}. (42) 

Thus, oriented hyperplanes P G are in a onc-to-onc correspondence 
with hyperplanes vr C M^'^ which make angle 7r/4 with the subspace 
= {(x,0)} CM^'^ 

(iii) An oriented hypersphere S C is in an oriented contact with an 
oriented hyperplane P C M"^, if and only if cr G tt. 

(iv) Two oriented hyperspheres 81,82 C are in an oriented contact, if 
and only if their representatives in the Minkowski space (Ti,(T2 G M^'^ 
differ by an isotropic vector: |(Ti — (T2I = 0. 

In the cyclographic model, the group of Laguerre transformations admits a 
beautiful description: 

Theorem 42 (Fundamental theorem of the Laguerre geometry) 

The group of Laguerre transformations is isomorphic to the group of affine 
transformations ofM.^'^: y 1— > XAy + b, where A G 0(A'", 1), A > 0, and 
6gM^'^ 
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